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Abstract. The abundance of local clusters is a traditional way to derive the amplitude of matter fluctuations, commonly speci- 
fied by <r 8 , but which suffers from a systematic uncertainty arising from the lack of accurate knowledge of the mass temperature 
relation. In the present work, by assuming that the observed baryon content of clusters is representative of the universe, we show 
that the mass temperature relation (M - T) can be specified for any cosmological model. WMAP constraints on the baryonic 
content of the Universe and the Sl M - H relation allows one further improvement in tightening this M - T relation. This 
approach allows to remove most of the above uncertainty, and to provide an estimation of cr 8 whose uncertainty is essentially 

■ statistical. The values we obtain are fortuitously almost independent of the matter density of the Universe (cr 8 ~ 0.6-0.63) with 
an accuracy better than 5%. Quite remarkably, the amplitude of matter fluctuations can be also tightly constrained to similar 
accuracy from existing CMB measurements alone, once the dark matter content is specified. However, the amplitude inferred 
in this way in a concordance model (A - CDM) is significantly larger than the value derived from the above method based 
on X-ray clusters. Such a discrepancy would almost disappear if the actual optical thickness of the Universe was but could 
also be alleviated from more exotic solutions: for instance the existence of a new non-baryonic light dark component in the 
Universe as massive neutrinos, with Sl d ~ 0.01 - 0.03. However, recent other indications of <r 8 favor a high normalization. In 
this case, the assumption that the baryonic content observed in clusters actually reflects the primordial value has to be relaxed 

either there exists a large baryonic dark component in the Universe with Q.J ~ 0.01 - 0.03 ~ 0.5f2i or baryons in clusters 
have undergone a large depletion during the formation of these structures. We concluded that the baryon fraction in clusters 
is not representative and therefore that an essential piece of the physics of baryons in clusters is missing in standard structure 
formation scenario. 
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1. Introduction measurements of clusters could in principle provide a direct 

measurement of this quantity but present-day results are con- 

The amplitude of matter fluctuations in the present-day uni- , ,. . rr, ■ .■ • ■ , 

r r ■> tradictory. This question remained unresolved because the am- 

verse is an important quantity of cosmological relevance. The ... , t .. . , c , , ... , , , .. 

r M J b plitude of matter obtained from clusters with hydrostatic equa- 

abundance of clusters is an efficient way to evaluate this quan- ^ leads tQ ^ v ^ ^ ^ Q J ± Q Q6 (Markevitch 1998> 

tity, commonly expressed by cr 8 , the r.m.s. amplitude of the ^ ^ ^ sdjak 2QQ2) M WMAp recently ofe _ 

matter fluctuations on the 8/z Mpc scale. A statistical preci- . ■ , A n , A , /c , . , OAAa s u ■ ■ , 

F F tamed cr g ~ 0.9 + 0.1 (Spergel et al. 2003). However, the vmal 

sion of a few % on cr« is possible from existing samples of , . ,. ffi . ■ , , • t * 

s F cluster masses are difficult to obtain, and values inferred with 

X-ray clusters, but in practice the relation between mass and j.^ . , £1 c ■ . 

J r different methods spread a large range of values; for instance 

temperature is needed for such evaluation: n , . , /OAAAN ■ . . t ^ , 

r Roussel et al. (2000) pointed out that hydrostatic mass estima- 
tions were lower than values inferred from numerical simula- 
tions and Henry (2004) recently found that published values of 
(Oukbir and Blanchard, 1992), Q M being the present-day mat- # a q uantit y proportional to A TU , could differ by a factor of 
ter density parameter and A being the contrast density relative nearl y two - This leaves a lar 8 e uncertainty on the actual ampli- 
to the Universe at the radius at which M 15 is taken. The value tude of matter fluctuations derived from clusters. In this paper, 
of At m has been estimated from X-ray properties of clusters by we P ro P ose a new a PProach to derive the mass-temperature re- 
different methods, essentially hydrostatic equations on one side latlon ln a self-consistent way. This allows us to combine the 
and numerical simulations on the other side, which lead to dif- bar y° n bud 8 et from the CMB ' observed gas mass in clusters 
ferent normalizations (from <r 8 ~ 0.6 to cr 8 ~ 1 .). Mass lensing and the P resent da y abundance of clusters to infer a tight con- 
straint on the amplitude of matter fluctuations obtained from 



T = A TM M 2 ' 3 (n M (l + A)/179) 1/3 /* 2/3 (l + z) keV (1) 
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Fig. 1. The central area is the amplitude of matter fluctuations 
expressed in term of <x 8 versus the normalization of the mass 
temperature Atm (Eq. Q for a flat model with H« = 0.3. 
Grey areas are our one, two and three sigma level contours. 
Filled squares are from Vauclair et al. (2003), open squares are 
Pierpaoli et al. (2001), around Atm ~ 7 and Pierpaoli et al. 
(2003), x symbol is from Evrard et al. (2002), Triangle is from 
Seljak (2002), inverted triangles are from Viana et al. (2002) 
and Viana et al. (2003). 
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Fig. 2. The red line is the central value of Atm which is the nor- 
malization of the mass temperature relation Eq.^ The WMAP 
relation between Hq and D.q has been used, as well as the con- 
straint on the baryon content of the Universe. One and two 
<r uncertainties arising from uncertainty on Q;, (= 0.023 + 
0.002) are shown as blue and yellow areas. Horizontal areas 
correspond to estimations of Atm from hydrostatic methods 
(light green) obtained by Roussel et al. (2000) and Markevitch 
(1998) and from numerical simulations (light blue) obtained 
from Bryan & Norman (1998) and Evrard, Metzler & Navarro 
(1996). 



the cluster abundance and compare it to the amplitude inferred 
from the CMB. 



2. Mass-Temperature relation 

2.1. The cr & -A TM degeneracy. 

The determination of <x 8 from the cluster abundance is a stan- 
dard procedure that has been used by many authors, leading to 
somewhat dispersed values. Here we use the Sheth and Tormen 
(1999) mass function and a sample of X-ray selected local clus- 
ters (f x < 2.210~' 1 erg/s/cm 2 and \b\ < 20 deg, Blanchard et al., 
2000 updated from BAX, Sadat et al. 2004). The relation be- 
tween <t%-Atm is presented in Figure [J for a flat model with 
Q.M - 0.3 with some other recent measurements, based on 
ROSAT samples of X-ray clusters and a recent analytical mass 
functions (del Popolo 2004). We do not include analyses based 
on HEAO-1 such as Henry (2004) or using the classical Press 
and Schechter expression as do Ikebe et al. (2001), or other 
mass functions. A NFW profile (Navarro et al. 1995) with c — 5 
was assumed when necessary. This shows that most of the dis- 
persion among different analyses (which use nearly the same 
clusters) is due to the different values used for the normaliza- 
tion constant A tm- Most of the remaining differences are due to 
differences in temperatures used (with or without cooling flow 
correction, or temperature cuts). Notice that the point present- 
ing the largest deviation (Viana et al. (2002)) is based on the 
luminosity function. We converted their mean luminosity to a 
mean temperature of 2.68 keV to derive an equivalent A TM . 



2.2. The baryon fraction argument 

Clusters are useful cosmological probes in several important 
ways. Their baryonic fraction ft can be inferred from observa- 
tions: 



fb = 



M,„, 



(2) 



the mass in observed baryons, M/,, consists mainly of the X- 
ray gas and of a small part of the stars (Roussel et al 2000), 
while the total mass M tot could be estimated through one of 
the above-mentioned methods. Under the assumption that the 
baryonic and dark matter amounts are representative of the uni- 
verse, the baryon fraction can be related to the cosmological 
parameters density Q./, and Q. m : 



f b = r 



(3) 



T is a numerical factor that has to be introduced in order to cor- 
rect for the depletion of gas during cluster formation and which 
can be determined only from numerical simulations (White et 
al. 1993). In practice, a good working value, at least in the outer 
part of clusters, is T = 0.925 (Frenk et al. 1999). The bary- 
onic content of the Universe is now known quite accurately 
through WMAP and other CMB measurements (u>b = Qbh 2 = 
0.023 + 0.002, Spergel et al. 2003; the statistical uncertainty 
being doubled in order to account for differences in various 
priors), essentially consistent with the abundance of Deuterium 
(Kirkman et al. 2003) and with the baryonic content of the IGM 
(Tytler et al 2004). 
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2.3. Self-consistent mass-temperature relation 

While the above relations have been widely used to obtain con- 
straints on Q. m assuming that the M - T relation is known, here 
we follow a different approach. Given the uncertainty in the 
actual value of A TM , we can use the knowledge of gas (and 
stars) masses in clusters and of Q/, to constrain the total mass 
in clusters, at least as a function of Q m and h, and thereby infer 
the mass-temperature relation. A slight source of complexity 
in gas mass measurements comes from the fact that the gas in 
clusters may be clumpy. If so, the gas mass estimation from 
average radial profile of the emissivity overestimates the actual 
gas mass by a factor C 1/2 , where C is a measurement of the 
clumping of the gas. Sadat and Blanchard (2001) have stud- 
ied in detail the change in shape of the gas fraction with radius 
in clusters: they found that the gas fraction follows rather well 
a scaling law, i.e. it is roughly identical among different clus- 
ters when expressed in term of the radius normalized to the 
virial radius. Furthermore they found that in the outer part the 
shape was close to what has been found in numerical simula- 
tions provided the outer amplitude is corrected for clumping 
(the value of T being roughtly constant ~ 0.925 for A < 1000). 
This implies that correction for clumping is indeed important 
to have an unbiased description of the internal structure of clus- 
ters. Mathiesen et al. (1999) found an average C 1/2 of 1.16 at 
the radius corresponding to a contrast density A of 500 (they 
also found that taking only clusters with no secondary peak at a 
level of 1 % of the global maximum reduced the average C 1 12 to 
1 .093). Because the clumping factor seems to vary rapidly with 
radius, it is safe to work on clusters at a similar radius. We have 
used the gas mass determination from Vikhlinin, Forman and 
Jones (1999, VFJ99 hereafter), using their most external radius 
for mass determination, which is /z-dependent, at the average 
temperature of 4 keV. VFJ99 provided gas mass measurements 
at the radius R1000 where the contrast density in the gas is 
1000 times some fiducial baryon density (2.85 10 9 M o /Mpc 3 ), 
which corresponds to nearly half of the best cjb derived from 
WMAP. Typical density contrasts at our working radius are in 
the range 480-625, at which we can directly apply the above 
correction for clumping. VFJ99 excluded clusters with double 
or very irregular X-ray morphology, a criteria that seems less 
demanding than the criteria for regularity used by Mathiesen 
et al. (1999). However, in both cases roughly one third of the 
clusters were excluded from the analysis. We therefore used 
the value C 1/2 = 1.093. The difference to C 1/2 = 1.16 is a 
source of systematic uncertainty on the final mass of 6%. We 
further corrected for a star contribution of 34h L5 % (Roussel 
et al. 2000). Knowing the baryon mass, relations 2 and 3 can 
be used to infer the total mass (depending on Q m ) at the radius 
R1000. Finally, in order to use the above mass estimation in the 
mass function we need to estimate the mass at the virial radius 
(unfortunately measurements of both apparent gas mass and 
clumping are not available at the virial radius). This virial mass 
can be estimated assuming a NFW profile with a fixed con- 
centration parameter c. Hereafter we used c = 5. WMAP has 
provided high precision data in a field where order of magni- 
tude estimations were the only possibility a few years ago. This 
allows us to constrain very tightly some quantities which are 



often combinations of a few cosmological parameters. For the 
present analysis, we use the location of the so-called Doppler 
peak which allows us to establish a tight relation between Q m 
and h in a flat universe (Page et al. 2003). The above procedure 
has been applied to derive the mass temperature normalization 
Aim as a function of Q m . The result is shown in Figure 1. As 
one can see the above procedure allows us to determine the 
value of Atm as a function of Q m with a small uncertainty : we 
found A TM ~ 4.9Q m 75 + 10%, for 0.3 < Q m < 1.0; the values 
cover the range of the various estimations based on the differ- 
ent approaches that we have mentioned. From this relation we 
can now infer the typical temperature of clusters formed from 
fluctuations within R = 8/T 1 Mpc spheres: 

7V>Mpc * A m (1.19Q m ) 2 /3 * 3.65keVQ m 009 (4) 

therefore, the amplitude of matter fluctuations cr 8 is essen- 
tially controlled by the abundance of clusters with tempera- 
tures around 3.5 keV almost independently of the value Q m 
(this comes from the fact that A TM varies with Q m and that 
this variation accidentally compensates almost exactly the vari- 
ation obtained for a fixed value of A TM )- Knowing the mass- 
temperature relation and its uncertainty we can determine the 
amplitude of matter fluctuations by fitting the local tempera- 
ture distribution function and assuming a F-like spectrum with 
F = 0.2 as explained in section 2.1. 

The result is shown in Figure 2a. As one can see, at a 
given value of Qm the amplitude of cr 8 is well constrained. 
Furthermore to the first order the best cr 8 is independent of 
Qm (o"8 ~ 0.63 ± 3.%(lcr) for Q A = 0.7. Interestingly this is 
close to the value obtained by Viana et al. (2002): cr 8 ~ 0.61). 
Our conclusion appears somewhat surprising as it differs from 
standard analyses based on a fixed normalization Atm, which 
cannot simultaneously account for the baryon fraction in a con- 
sistent way for arbitrary Q M . Gas masses from the VFJ99 sam- 
ple present a moderate dispersion of the order of 20% (Sadat et 
al. 2005), implying rather small uncertainties on our gas frac- 
tion estimates of the order of about 4% at our working radii, 
which will produce an uncertainty on cr 8 of 2.5%. More impor- 
tant is the correction for clumping. For instance, Voevodkin & 
Vikhlinin (2004) have estimated cr 8 from the baryon mass func- 
tion in a Cold Dark Matter framework. In the case Q« ~ 0.3 
their approach is very close to ours, but they used gas mass es- 
timation at the virial radius and did not correct for clumping. 
This leads to virial masses which are ~ 20% lower than ours, 
leading to A TM ~ 1 1 keV and therefore cr 8 ~ 0.7, in very good 
agreement with their estimation. 

3. Need for Dark Matter 

The amplitude of matter fluctuations is strongly constrained by 
the CMB data. In the following we use the constraint on cr 8 in 
a concordance model obtained from the CMB fluctuation anal- 
ysis including the temperature-polarization cross power spec- 
trum (TE) by the WMAP team (Kogut et al 2003). 

The comparison of the value of <x 8 from CMB data with 
the one from clusters reveales a critical discrepancy between 
the two measurements (Figure 2a). It is clear that within any 
model with Q A ~ 0.7 the amplitude of <x 8 we derived from 
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Fig. 3. a) top: The amplitude of matter fluctuations from clus- 
ters abundance using the mass-temperature relation found in 
the present analysis compared to the amplitude of matter fluc- 
tuations derived from CMB data (Grainge et al 2003, Pearson et 
al. 2003, Ruhl et al. 2003, Kuo et al, 2004). The grey area corre- 
sponds to 1, 2, 3 cr contours on two parameters, dashed lines are 
contours on one parameter. The one and two sigma amplitudes 
obtained from an average of recent weak shear measurements 
are also shown as dashed regions (see text for references), b) 
bottom: amplitude of matter fluctuations from CMB data ver- 
sus optical thickness. 



clusters, cr% = 0.63 + 0.02, is significantly smaller than what 
is expected from the CMB alone (cr 8 = 0.88 + 0.035), which 
is close to the recent determination from the combination of 
WMAP and the Lyman a forest (Seljak et al 2004). 

The non-zero optical thickness t, which is requested only 
from the TE spectrum, is a key factor in this discrepancy: the 



high erg obtained from CMB data depends critically on the ac- 
tual amplitude of the optical thickness r (see Figure 2b), al- 
though forcing r = does not entirely remove the discrep- 
ancy. An accurate knowledge of t is therefore critical to prop- 
erly evaluate the amplitude of matter fluctuations in the concor- 
dance model. One can see from Figure 2b how much the value 
of erg obtained from CMB data depends on the actual value 
of the optical thickness and remains the main source of uncer- 
tainty in establishing the value of cr 8 . We have also checked that 
when CMB data are restricted to the range 400 < I < 1200, the 
above discrepancy remains essentially unchanged. Allowing a 
non-power law initial power spectrum is therefore not expected 
to solve this issue. 

Here above, we have considered models in which the dark 
matter is only made of cold dark matter, the dark energy be- 
ing a pure cosmological constant (in terms of the equation of 
state of vacuum p = wp, this means w = -1), and that X-ray 
gas and known stars are the only existing baryons in clusters. A 
first possibility to investigate is to examine whether a different 
equation of state for the vacuum, so-called quintessence, might 
solve this discrepancy. We have therefore investigated flat mod- 
els with arbitrary w and quintessence content Q.q. Indeed com- 
binations of CMB and cluster data are known to provide tight 
constraints on such models (Douspis et al. 2003). With the ap- 
proach developed here, models which were found to match 
CMB and clusters were found to satisfy the following con- 
straints: 0.46 < Q.q < 0.54 and -0.5 < w < -0.4. Such 
models are currently at odds with constraints on quintessen- 
tial models (Douspis et al. 2003; Tegmark et al. 2003; Riess 
et al. 2004) resulting from the combination of various data in- 
cluding type la supernovae data. We therefore require an al- 
ternative approach to solve the above issue. In the following, 
we examine whether the introduction of an additional compo- 
nent of the dark matter content of the universe would remove 
the above discrepancy. Neutrinos are known to exist and to be 
massive, so perhaps the most natural massive component of the 
universe to be introduced is in the form of a neutrino contribu- 
tion. This solution has already been advocated to solve this dis- 
crepancy in an Einstein de Sitter Universe (Elgar0y & Lahav 
2003, Blanchard et al. 2003). Indeed, the presence of a light, 
but non-zero, component of dark matter significantly modifies 
the transfer function of primordial fluctuations which results 
in a lower amplitude on small scales. Given existing measure- 
ments of mass differences we consider only the case where 
the masses are equal. Within a concordance model (Qa = 0.7 
Q. m = 0.3), by combining the constraints from CMB and clus- 
ter data, and marginalizing on (u>b, Hq, n, t) we found that a 
contribution of Q v = 0.016 ± 0.003 is preferred with a sig- 
nificance level well above 3<x (see Figure 3a), improving the 
significance of such possible evidence compared to Allen et 
al. (2003). This confirms that the presence of a small contribu- 
tion of neutrinos, with a typical mass of .25 eV, to the density 
of the universe allows one to reconcile the amplitude of mat- 
ter fluctuations from clusters with the one inferred from CMB 
data. We notice that such a value is above the upper limit in- 
ferred by the WMAP team using a combination of several as- 
tronomical data (Spergel et al. 2003). Finally, weak shear esti- 
mations have provided measurements of the amplitude of mat- 
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ter fluctuations which can be compared to that obtained from 
clusters (Refregier 2003). There are some differences in pub- 
lished values which probably reflect systematic uncertainties 
not yet fully identified. However, taking the independent mea- 
surements of erg from weak lensing obtained from an average 
of recent measurements (Bacon et al. 2003, Brown et al. 2004, 
Chang et al. 2004, Hamana et al. 2003, Heymans et al. 2004, 
Hoekstra et al. 2002, Jarvis et al. 2003, Massey et al., 2004, 
Refregier et al. 2002, Rhodes et al. 2004, Van Waerbeke et al., 
2004) which lead to an acceptable^ 2 , from WMAP and Lyman 
a forest (Seljak et al 2004), and the value of /3 from 2dFGRS 
(Hawkins et al. 2003), the low amplitude of cr 8 obtained above 
is not favored. We are therefore left with the conclusion that our 
initial assumption that baryons in clusters are fairly representa- 
tive of baryons in the universe is unlikely and therefore that the 
observed amount of baryons in clusters does not reflect the ac- 
tual primordial value (a possibility that has been advocated by 
Ettori 2003). Several mechanisms could lead to this situation: 
the most direct way could be the fact that a significant fraction 
of the baryons are in a dark form, either in the Universe or in 
clusters (for instance either in the form of Machos, or in a large 
gaseous unidentified component, Bonamente et al, 2003), or 
that a significant fraction of the baryons has been expelled from 
clusters during their formation process. In such cases, the ob- 
served Mi, is biased to lower values. The actual mass of clusters 
from Eq. 2 and Eq. 3 can then be obtained assuming a depletion 
factor 1 - / implying that /O/, represents the missing baryons. 
Again the combination of CMB and cluster constraints allows 
us to evaluate the amplitude of /O/,. From Figure 3b, one can 
see that such a component, /Oj ~ 0.023, should represent 
nearly half (T ~ 0.5) of the primordial baryons in order to 
solve the discrepancy. Although heating processes are required 
to account for the observed properties of X-ray clusters, they 
currently do not lead to such a high level of depletion (Bialek, 
Evrard&Mohr2001). 

4. Conclusions 

The determination of the amplitude of matter fluctuations 
within pure cold dark matter model, using two methods, 
namely the CMB and the local cluster abundance, leads to two 
significantly different values. There are several ways to solve 
this discrepancy, although each represents a noticeable depar- 
ture from the standard concordance model. The existence of 
a non-baryonic dark component, like a neutrino contribution, 
would allow us to solve this discrepancy, although such a so- 
lution leads to a low value of <x 8 which is not favored by other 
evidence. If the actual value is larger, <x 8 ~ 0.8 - 0.9, the un- 
avoidable conclusion is that the baryonic content of clusters at 
A ~ 500 is not representative of the Universe. In this case, an 
astrophysical solution could be that baryons in clusters could be 
in a dark form, or at least undetected until now. Alternatively, 
baryons in clusters could have been severely depleted implying 
that the actual value T is much smaller than the value we used 
above, the apparent baryon fraction being biased to low values 
compared to the actual primordial value. Finally, several obser- 
vations might help to clarify this issue: the above conclusion re- 
lies on the actual value of the optical depth r found by WMAR 
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Fig. 4. Constraints on O v and Qdarkbarytmh 2 given by the com- 
bined analysis of CMB and Cluster data. The dark energy com- 
ponent Qa has been set to 0.7 in a flat cosmology. 



If the actual value was consistent with zero most of the discrep- 
ancy would disappear. Confirmation of the actual value of t is 
therefore critical and its better estimation will allow a better 
estimation of cr 8 from the CMB. Other sources of information 
on o"8 will also obviously clarify this issue: weak lensing can 
potentially allow one to directly measure the actual amplitude 
of matter fluctuations with a similar precision to what has been 
obtained here with clusters, provided that systematic uncertain- 
ties are fully understood; the clusters masses could be mea- 
sured from their lensing signal providing a direct estimation 
of the normalization constant A TM , allowing one to distinguish 
between low and high normalizations. Other direct measure- 
ments of the amplitude of matter fluctuations like those derived 
from the Lyman-ff forest power spectrum (Croft et al. 1998) 
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could also help clarify this issue. It is remarkable that some of 
the observations that are expected in the near future can po- 
tentially bring fundamental information on clusters physics or 
alternatively may reveal the existence of a previously unidenti- 
fied type of dark matter with Q^m as low as 0.01. 
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